The coordination dependence of electric-field gradients and hyperfine fields 
at 5sp impurities on fee metal surfaces 
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We present a comparison between accurate ab initio calculations and a high-quality experimental 
data set (1986-2002) of electric-field gradients and magnetic hyperfine fields of Cd at different sites 
on Ni, Cu, Pd and Ag surfaces. Experiments found a systematic rule to assign surface sites on 
(100) and (111) surfaces based on the main component of the electric-field gradient, a rule that does 
not work for (110) surfaces. Our calculations show that this particular rule is a manifestation of 
a more general underlying systematic behavior. When looked upon from this point of view, (100), 
(111) and (110) surfaces behave in precisely the same way. The experimentally observed parabolic 
coordination number dependence of the Cd magnetic hyperfine field at Ni surfaces is verified as a 
general trend, but we demonstrate that individual cases can significantly deviate from it. It is shown 
that the hyperfine fields of other 5sp impurities at Ni surfaces have their own, typical coordination 
number dependence. A microscopic explanation for the different dependencies is given in terms of 
the details of the s-DOS near the Fermi level. 

PACS numbers: 68.47.De, 68.35. Fx, 68.35.Dv, 75.70.Rf, 76.80.+y, 71.20.-b 



I. INTRODUCTION 



Nuclear probe techniques, such as Mossbaucr spec- 
troscopy, Perturbed Angular Correlation spectroscopy 
(PAC), Nuclear Magnetic Resonance (NMR) and oth- 
ers, have proven in the past to be versatile tools for the 
study of a broad set of phenomena. By simultaneously 
detecting both electric-field gradients (EFG) and mag- 
netic hyperfine fields (HFF), it has been possible to ob- 
tain a finger print of the electronic and magnetic config- 
uration near and at the probe nucleus Ab-initio 
band structure calculations have assumed an important 
role in supporting, improving or confuting the interpre- 
tation of experimental data. Especially for such sensitive 
properties as HFF's or EFG's, the availability of these 
parameter-free calculations with an accuracy that is com- 
parable with the experimental error bars, has opened new 
possibilities. For instance, the general behavior of the 
HFF induced on sp- and d-impurities embedded in bulk 
magnetic materials, especially in Fc and Ni, has been well 
characterized and the calculated values compared fairly 
well with the experimental onesiSiSiiSili Similarly, the 
systematics of the EFG on 5sp-impurities in non-cubic 
metals (Zn, Cd, Sb) have been studied by experiments 
and calculations (Ref. and references therein). Se- 
lected calculations have been performed on more exotic 
systems. To name a few, the HFF and EFG on a Cd 
impurity atom at Fe/Coi^and Fe/Agi^i^ interfaces have 
been studied, focussing on the relation between the HFF 
induced on the radioactive probe atom and the magnetic 
moment profile across the interface. Similar studies have 
been performed on magnetic nanoclusters embedded in 
Agi& and Cui£ matrices. In thin layers of fee Fe on a 
Cu substrate, HFF and EFG at the surface and interface 



have been calculated^ 

If, for a few host materials, the dependence of HFF 
or EFG on the atomic number of the impurity atom has 
been studied in detail, less investigated is the variation in 
these properties for one specific probe atom in different 
environments derived from the same host material. An 
example of such a study - motivated by the unique capa- 
bility of nuclear probe techniques to follow the diffusion 
of single atoms - is provided by the experiments of Voigt 
et al.j£ recently extended and completed by Potzger et 
al4£ Those authors put a Cd atom at different low-index 
Ni-surfaces and at kinks and steps on those surfaces, and 
measured the HFF and EFG at Cd. Overlooking the now 
fairly complete data set, Potzger et al. conclude that the 
HFF at Cd depends more or less parabolically on the 
number of Ni-atoms in the first nearest neighbor shell 
(NN-shell) (black dots in Fig. |3J), and not on the exact 
position of those Ni-neighbors (dubbed 'symmetry inde- 
pendence of the HFF' in Ref. l20f) . This is attributed to 
a gradual change in the local DOS, not further specified. 
One of the two main goals of the present paper is to assess 
the validity of this claim, to elucidate the possible physi- 
cal mechanism behind it (Sec. lIV"H) . and to try to extend 
this 'rule' to impurity atoms other than Cd (Sec. IIV"2|l . 
all this by calculating the HFF and EFG at Cd in dif- 
ferent Ni-environments by ab initio methods. A second 
goal is to examine the generality of common guidelines 
for site identification of Cd on (100) and (111) metal 
surfaces based on EFG measurements, as were estab- 
lished from experiments i^^i We will demonstrate that 
these guidelines are a manifestation of a simple and gen- 
eral physical mechanism that is valid for all low-index 
surfaces (Sec. iHljl . Our work follows pioneering cluster 
calculations by B. Lindgren et ali 2 ^ 2 ^ 2 ^ 2 ^, who investi- 
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gated HFF and EFG of Cd as adatom or in a terrace 
site at (100) and (111) Ni surfaces, and a more system- 
atic work by Ph. Mavropoulos et al3i for probe atoms 
belonging to the whole 4sp series (Cu to Sr), placed on 
Ni and Fe (100) surfaces. Independent from our work, 
some of the questions that will be discussed here were 
studied very recently by Ph. Mavropoulos^l using a dif- 
ferent method (full-potential KKR) , a different exchange- 
correlation functional (LDA) and not considering atomic 
relaxations. 



II. METHOD AND DETAILS OF THE 
CALCULATIONS 

For the present calculations we have employed state of 
the art first-principles techniques, developed within the 
Density Functional Theory (DFT)2S*2a. Most of the cal- 
culations have been performed using the full-potential 
augmented plane wave + local orbitals (APW+lo) 
method as implemented in the WIEN2k packaged We 
have simulated the surface+impurity complexes by using 
the slab-supercell approach. The impurities have been 
put both in an adatom site and at a substitutional po- 
sition in the surface layer (called the terrace site from 
now on) at the three low-Miller-index surfaces, which 
in turn were composed of repeated slabs of 5 to 7 sub- 
strate layers (up to 11 for test calculations), separated 
by some vacuum space. The size and shape of the cells 
were chosen such to avoid artificial interaction between 
the impurity atoms. For symmetry reasons, the impu- 
rity atoms were placed at both sides of the slabs. Test 
calculations showed that about 6 A of vacuum (equiva- 
lent to about 3.5 atomic layers) is sufficient to decouple 
the Cd atoms at the two surfaces, and that a similar dis- 
tance is needed to decouple the Ni atoms within the same 
plane. More details of the slabs and the vacuum are given 
in tables later in the paper. As exchange-correlation 
functional, the PBE generalized gradient approximation 
(GGA)2i was used. A few test calculations (see later) 
have been performed with the local density approxima- 
tion (LDA). In the APW+lo method, wave functions, 
charge density, and potential are expanded in spherical 
harmonics within nonoverlapping atomic spheres of ra- 
dius Rmt and in plane waves in the remaining space of 
the unit cell. Rmt values of 2.2 a.u. for Ni and 2.4 a.u. 
for the impurity were chosen throughout. The maximum 
I for the waves inside the atomic spheres was confined 
to i m ax = 10. For the majority of the cases, the wave 
functions in the interstitial region were expanded in plane 
waves with a cutoff of K mo:c =7.0/R™y = 3.18 a.u. -1 or 
better, while the charge density was Fourier expanded 
up to G mQX =16 VRy. For a few large cells, a lower 
Kmaa;=2.73 a.u. -1 was taken, which was tested to yield 
still reasonable forces. For the Brillouin zone sampling, 
a mesh of special k points corresponding to at least a 
14 x 14 x 14 mesh (104 points) for bulk fee Ni was taken. 
The lattice constant of Ni used throughout the calcu- 



lations is the one which minimizes the GGA bulk total 
energy, i.e. a%^ A = 3.51 A. For Cu, Pd and Ag, the 
experimental lattice constants were used. 



When the perfect crystal structure is broken, as in 
presence of an impurity, atomic relaxations of the atoms 
in the vicinity of impurity appear. Such relaxations have 
rather long range character, and lead to a volume change 
of the crystal, as testified by refined lattice parameter 
measurements. Additionally, since they alter the nature 
and strength of bonds, they change local properties such 
as HFF and EFG as well. Atomic relaxations tend to be 
the larger in size the more open is the crystal structure 
(e.g. zincblcnd) 3 ^, but it has been shown in the litera- 
ture to be important also in close-packed matrices for 
HFF'o 10 ! 11 and even more for EFG's^ 3 .. As taking into 
account atomic relaxations is computationally quite ex- 
pensive, we have speeded up the calculations by using 
a combination of methods. First, the atomic positions 
were relaxed using the pscudopotential plane wave VASP 
cod o 34 ' 35 . The all-electron APW+lo code was then used 
in a second stage to further relax the atoms to their equi- 
librium positions (with forces less than 1 mRy/a.u.), and 
to calculate the hyperfme properties. 



Before proceeding with the discussion of the results, 
it is important to assess the accuracy of our calculated 
HFF's or EFG's. To this end, let us look at a well-known 
case, Cd in bulk Ni. The experimental value is -6.9 T (4.2 
K) 3 ^. Our calculation, carried out for a large 2x2x2 
32-atom supercell with relaxed nearest neighbors, yields 
-12.0 T. The calculation included SO-coupling (in con- 
trast to all other calculations reported in this paper), 
from which we could obtain the orbital and dipolar con- 
tributions to the hyperfme field. They were smaller than 
0.01 T, however, as could be expected for a closed shell 
atom as Cd ([Kr]5s 2 4d 10 ), such that the so-called Fermi 
contact term^I is the only contribution to the HFF. The 
Ni-moment was put in the [111] direction, as in nature. 
Without relaxation, the HFF was -9.8 T. Having taken 
all these precautions, we must attribute the 5.1 T dif- 
ference with experiment entirely to the approximations 
contained in GGA. This magnitude of error is similar to 
the 4.2 T overestimation of the HFF for Cd in Fe, and 
is of the same magnitude as was found for all 5 th period 
impurities in Fe>i2*ii This is the state-of-the-art for the 
accuracy that can be obtained currently for calculated 
hypcrfinc fields. In what follows, we can therefore ex- 
pect calculated hypcrfinc fields that are 5 T too negative 
if this error is independent on the coordination number, 
or alternatively - if the error is more random and does 
not depend on the coordination number - we have to put 
error bars of 5 T on our results. Electric-field gradients 
instead can be calculated with more accuracy, and we 
expect here an error bar better than (±2 ■ 10 21 V/m 2 ) 
(see for instance Ref . l38j) . 
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III. ELECTRIC-FIELD GRADIENTS 

For an efficient discussion later in this section, we first 
summarize a few properties and conventions with respect 
to the electric-field gradient (EFG) at a nucleus in a 
crystal (for a more elaborate discussion, see for instance 
Refs. and liflh . The EFG is a traceless, symmetric 
tensor of rank 2 (5 components), and is a measure of 
the deviation from spherical symmetry of the electronic 
charge density around the nucleus under consideration. 
In its principal axis system (PAS), the EFG tensor is 
diagonal with diagonal elements V xx , V yy and V zz and 
V KX +V yy +V zz = (traceless, hence specifying 2 out of the 
3 quantities is enough) . The axes of the PAS are conven- 
tionally labelled such that \V ZZ \ > \V yy \ > |Vx X |. Instead 
of specifying 2 of the 3 diagonal elements, usually the 
couple (V zz , rj) is specified, with the so-called asymme- 
try parameter which equals (Vxx — V yy )/V zz (0 < r] < 1). 
Often one refers with the term 'electric-field gradient' to 
V zz only. Whenever we use the word 'EFG', we mean 
the tensor with its 5 components. If we mean the main 
component in the PAS, we write V zz . 

Over the past two decades, a large body of experi- 
mental results has been obtained for the EFG's of Cd 
at different surfaces of Ni, Cu, Pd and Ag (see Tab. U 
for a summary of the experimental work, with references 
tracked down as much as possible to the original sources). 
In contrast to what the structure of Tab.|l]might suggest, 
the assignment of the experimental values to either ter- 
race or adatom sites does not happen in a direct way. For 
(100) and (111) surfaces, it is easy to separate terrace and 
adatom sites from all other ones (steps, kinks) due to the 
axial symmetry {v = 0) of the EFG and the perpendicu- 
lar orientation of the z-axis of its PAS. It is much harder 
to deduce which of both is the terrace site and which one 
is the adatom site. Only for a few cases a thorough ex- 
perimental analysis has been performed: Ag(100)2i and 
Pd(lll)22i^. Arguments include the thermal annealing 
behavior, the fraction of different sites dependent on the 
type of vicinal cut, and not at least a comparison with 
early ab initio cluster calculations Apart from these 
two well-tested cases, the partial information from ther- 
mal annealing behavior only (Cu(lll)i^) or from agree- 
ment with cluster calculations only (Cu(100)^) further 
supports the assignment in Tab. [I] To our knowledge, 
the results on Pd(110) are published only recently and 
in an indirect way (unpublished Ref. l46l mentioned in 
Ref . l25l) . Prior to this information, Tab. [I] was commonly 
summarized by the following practical 'rule': on a metal 
surface, a large V zz for Cd indicates the terrace site, while 
a low V zz indicates the adatom site. This rule has been 
applied e.g. by Potzger et ali^ to assign the adatom and 
terrace sites on Ni(100) and Ni(lll). If the rule would be 
not universal, this site assignment could be wrong, and 
consequently also the conclusion on the parabolic coor- 
dination number dependence of the Cd-HFF (Fig. [3] and 
Section llV lfl . Therefore, we will first examine the valid- 
ity of this rule. Is it really valid for all (100) and (111) 



TABLE I: Experimental absolute values of the principal com- 
ponent |14z| of the electric-field gradient tensor of Cd at var- 
ious metallic surfaces, with references tracked down as much 
as possible to the original sources. Units: 10 21 V/m 2 . The 
lowest available measurement temperatures were chosen, and 
are specified below (0 K if a reliable extrapolation was avail- 
able) . 
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"unpublished Ref. l4ll mentioned in Refs. Eland 13 
~_T=270-300 K 

T=0 K 
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I! T=0 K 
^unpublished Ref. |41| 
9 unpublished Ref. |41] 
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mentioned in Ref. V = 0.97 
— ia ,, - „. u>4f =77 K 
"unpublished Ref. |4fl mentioned in Ref. |23 r\ = 0.42 
°Ref.ELT=0 K 
PRef.|ll T=34QK 
«RefsT47j an d lip. T=0 K 
r Ref.p and Hj T=77 K 
s Ref. 4c , T=0 K 
'Ref. It , T=36 K 
"Ref. El and EJ T=77 K 



surfaces of the nonmagnetic metals Cu, Pd and Ag (i.e. 
also for the cases that are not experimentally tested in 
detail)? Can the rule be transferred without changes to 
the magnetic metal Ni? And is it really true only for 
(100) and (111) surfaces, and not for (110) (as the recent 
Pd information suggests)? Finally, what is the physical 
mechanism that lies behind this rule? 

In response to these questions, we present in Tab. HIT a) 
the calculated V zz for both relaxed and unrelaxcd struc- 
tures, together with their assigned experimental values 
(taken from Tab.[Q). We did not find in the literature any 
experimental investigation for Cd at the Ni(110) surface. 
Instead, for a qualitative comparison, the data for Cd at 
the Pd(110) surface are enclosed in parenthesis in Tab. 
IrETah We also include the results obtained with the clus- 
ter method by B. LindgreniS Unfortunately, the cluster 
calculations of Lindgrcn spanned all the orientations only 
for the Cu surface, while for Cd on Ni only the terrace 
sites at Ni(100) and (111) surfaces have been simulated. 
Results for Cu and Pd surfaces are therefore included in 
Tab. [nia) in parentheses as well. We recall here that 
differently to the cluster results of Lindgren, our calcu- 
lations arc for bulk slabs and not for non-periodic clus- 
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TABLE II: (a) Column 1: unrelaxed slab calculations for Vz Z 
of Cd on Ni surfaces. Column 2: idem, but relaxed. Col- 
umn 3: unrelaxed cluster calculations for V Z z of Cd on Ni 
surfaces, when available. If no values were available, results 
for Cu or Pd surfaces are given instead. Column 4: Abso- 
lute value of the experimental V Z z for Cd on Ni surfaces, if 
available (if not, values for a Pd surface are given), (b) Re- 
laxed slab calculations for Kz of Cd in terrace and adatom 
sites at one magnetic and three nonmagnetic (100) surfaces 
(sign included) , compared with the absolute value of the cor- 
responding experimental V zz . Units: 10 21 V/m 2 everywhere. 



(a) 



APW+lo 


Cluster Calc. 


|Exp.| 


Unrelax. 


Relax. 


Ref.[25j 


(see Tab. GJ 


(100) Terrace 10.5 


8.6 


11.0 


8.2 


(100) Adatom 4.0 


1.9 


(-0.8)cu 


2.8/0.3 


(110) Terrace 9.5 


-9.7 


(11.2) Pd 


(7.9) Pd 


(110) Adatom -9.3 


-6.4 


(-8.3) Pd 


(8.5) Pd 


(111) Terrace 10.8 


13.0 


16.5 


11.5/12.3 


(111) Adatom -1.4 


-2.9 


(-6.2)cu 


1.0 


(110) Terrace rj = 0.83 


rj = 0.74 


7?Pd = 0.97 


?7Pd = 0.97 


(110) Adatom rj = 0.18 


rj = 0.23 


77 Pd = 0.28 


J?p d = 0.42 


(b) 


(100) terrace 


(100) adatom 


APW+lo 


|Exp.| 


APW+lo 


|Exp.| 


Ni 8.6 


8.2 


1.9 


2.8/0.3 


Cu 9.9 


10.3 


-2.0 


0.8 


Pd 9.1 


8.2 


3.4 


2.8 


Ag 8.5 


7.5 


0.4 


0.3 



ters. Moreover we do a 'full' relaxation of all the atomic 
positions in the cell, while due to insufficient accuracy 
of the method, no real forces have been calculated by 
Lindgrcn: the distance between Cd and the surface was 
fixed at the experimental distance, and the surrounding 
Ni atoms were kept in their ideal fee positions (in con- 
trast to this, the Cd atom in our unrelaxed calculations 
was at an ideal fee lattice site). Our fully relaxed results 
compare generally very well with the experiments, pre- 
dicting for (100) and (111) a low V zz (close to zero) for all 
adatom sites, and a considerably higher V zz for the ter- 
race site, in compliance with the experimental rule. For 
the (110) surface, terrace and adatom sites have a com- 
parable V zz , again in agreement with the experimental 
observation. Taking into account the rather different ap- 
proximations, there is also reasonable agreement between 
our unrelaxed slab calculations and the unrelaxed cluster 
results, the former being generally closer to experiment 
(especially for (111)). For the (HO)-surface, the EFG is 
not axially symmetric (rj =/= 0). Values for 77 are included 



in Tab. [H{a) as well. Relaxation can be seen to have a 
moderate influence on rj, and the considerably larger rj in 
terrace sites as seen in experiments, is well-reproduced. 
Direct comparison with experimental values for rj is diffi- 
cult with the present data set, however (no experiments 
for Ni(110) and no relaxed calculations for Pd(110)). 

For the (110) terrace site we observe a change in the 
sign of V zz , when moving from the unrelaxed to the re- 
laxed system. This is readily explained as follows. V zz 
is defined as the largest component of the electric field 
gradient tensor in its principal axis system. But it might 
happen that there are two components of similar size, 
that have then necessarily opposite signs and yield a large 
asymmetry parameter rj (=close to 1). This is exactly 
what we found for the (110) terrace case. There is a 
large positive component along the normal on the (110) 
surface, while an almost as large negative component ex- 
ists along the [110] direction in the (110) plane. Upon 
relaxation the latter negative component becomes in ab- 
solute value larger than the positive one, and the latter 
axis becomes the one along which V zz is defined. In the 
adatom case, the z-axis is the [001] direction in the (110) 
plane. Experiments on the Ni(110) surface predict in- 
stead the EFG z-axis to lay along the normal to the (110) 
plane. This apparent deviation between experiment and 
calculation is only a consequence from the way how V zz 
is defined, not from a real discrepancy. For the other two 
substrate orientations, (110) and (111), the asymmetry 
parameter rj is zero and there is therefore no doubt on 
which component is the largest one: we find in all cases 
the principal axis along the surface normal, as in experi- 
ment. 

Having verified that the practical site assignment rule 
is valid for magnetic Ni(100) and Ni(lll) surfaces, and 
not for Ni(110), we now examine for one type of surface 
- (100) - whether other elemental metal substrates show 
the same behavior. Tab. ITlTb) shows the calculated and 
experimental V zz at the (100) surfaces of Ni, Cu, Pd and 
Ag. The calculations are in very satisfying agreement 
with the experiments, and it is safe now to conclude that 
the EFG in all of these 24 surface sites behaves in the 
same way. 

Finally, we search for the physical mechanism that is 
responsible for this systematic behavior. Instead of con- 
sidering for this the Ni-calculations from Tab. [HI we cal- 
culated the EFG of Cd at a complete set of unrelaxed Pd- 
surfaces. This avoids the complication of an EFG coming 
from spin up and spin down electrons, and it allows for a 
more direct comparison with experiment because for Pd 
all six cases are measured (Tab.^J). It does not harm not 
to consider relaxation here, because now we are looking 
for a global mechanism, not for fine details. The results 
are given for every diagonal component of the EFG sepa- 
rately in Tab. Hill together with the corresponding exper- 
imental components (mind the unconventional labelling 
of the axes, as specified in the caption of Tab. lIIIfl . Con- 
sidering the fact that these arc unrelaxed calculations, 
they reflect the experimental trend very well and hence 
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TABLE III: The three calculated (Vj; alc , unrelaxed) and experimental (V^? xp , see Tab.QJfor references) components of the EFG 
tensor in its PAS for Cd at 6 Pd surface sites. Units: 10 21 V/m 2 . The usual convention for labelling the axes of the PAS is not 
followed: the z-axis is perpendicular to the surface, x- and y-axes are in the surface plane. For (100) and (111) surfaces, all 
directions in the surface plane are equivalent for the EFG. For the (110) surface, our x-axis refers to the [llO] direction, and 
our y-axis to the [001] direction. The principal component - which would be labelled Vzz in the usual convention - is printed 
in bold. The sign of the experimental values was not measured, but is chosen here to agree with the calculations (for the small 
numbers for NN=3, this is nothing more than a guess). The number of Cd-5p electrons split into np alc , npy lc and np alc is given 
as well. 





(t NN 


t ycalc 


v yy 


t 7"calc 
* zz 




T/oxp 

v yy 
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v zz 
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(Ill) terrace 


9 


-4.8 


-4.8 


9.6 


-5.1 


-5.1 


10.2 


0.1478 


0.1478 


0.0952 


(100) terrace 


8 


-4.4 


-4.4 


8.8 


-4.1 


-4.1 


8.2 


0.1321 


0.1321 


0.0876 


(110) terrace 


7 


-7.6 


-0.5 


8.1 


-7.8 


-0.1 


7.9 


0.1214 


0.0986 


0.0736 


(110) adatom 


5 


2.0 


-6.2 


4.2 


-2.5 


-6.0 


8.5 


0.0681 


0.0972 


0.0651 


(100) adatom 


4 


-1.4 


-1.4 


2.8 


-1.4 


-1.4 


2.8 


0.0645 


0.0645 


0.0580 


(111) adatom 


3 


0.4 


0.4 


-0.8 


0.2 


0.2 


-0.4 


0.0424 


0.0424 


0.0574 



these Pd-surfaces form a reliable model system to search 
for the physical mechanism. 

Already in the very first successful ah initio calcula- 
tions of EFG's in metals, it has been shown that mag- 
nitude and sign of the principal component V zz are re- 
flected in the so-called asymmetry count An p of the p- 
electrons (for Cd, the completely filled d-shell does not 
contribute to the EFG). The asymmetry count is defined 
as An p = \ (n Px + n Py ) — n p „, with n Pi the number of 
electrons in the p ,;-orbital (see Ref. for the first analy- 
sis, and e.g. Refs. I50l5ll52l53| for applications): spherical 
symmetry (n Px = n Py = n Pz ^ leads to An p = (hence 
V zz = 0), charge accumulation along the z-axis (n Px 
large) leads to An p < (hence V zz < 0), etc. For our 
analysis, we slightly extend this idea and take as a work- 
ing hypothesis that sign and magnitude of any diagonal 
component Va is reflected in a generalized asymmetry 
count An l p = ^ (n Pj + n Pk ) — n Pi (and cyclic permuta- 
tions). This hypothesis is checked by Fig. ^a), where 
the calculated Va from Tab. Unl (right scale) are com- 
pared with the asymmetry counts An p derived from the 
calculated number of pi electrons in Tab. IIIII (left scale). 
With the right choice of scales in Fig. D^a) , the An p and 
Vii data points almost coincide, which means that apart 
from a unique scaling factor both quantities reflect ba- 
sically the same behavior. This proves our working hy- 
pothesis: the observed behavior of the components of the 
EFG is nothing more than a manifestation of the prop- 
erties of the distribution of the Cd 5p-electrons over the 
x-, y- and z-directions. Hence, if we understand the be- 
havior of the n Pi in Tab. IIIII we understand the EFG. 
This greatly facilitates our search, as - in contrast to the 
Va - the n Pi have a direct and intuitively understandable 
relation to the chemical bonds between Cd and its Pd 
neighbors. 

The coordination dependence of the n Pi is visualized in 
Fig. ^b), together with a least squares fit. The picture 
is complicated a bit by the fact that n Px and n Py are not 
identical to each other for (110) surfaces. To make the 



discussion more transparent, we will therefore first as- 
sume a hypothetical (110) surface with axial symmetry 
(?/ = 0) for which n Px = n p (gray symbols in Fig.^ they 
are the average of the calculated n Px and n Py ). We now 
clearly sec an almost linear dependence of n Px (or n Py ) 
and n Pz on NN (least squares fit in Fig.^b)): both n Px 
[n p ) and n Pz decrease linearly while reducing the coordi- 
nation, the decrease being stronger for n Px (n Py ) than for 
n Pz . If we calculate the An l p corresponding to the least 
squares fits in Fig. ^b) , we obtain the straight lines in 
Fig. ^a) (for overall axial symmetry, gray symbols for 
NN=(5,7)): this proves that the two-slope model from 
Fig. ^b) is indeed the key responsible for the observed 
behavior of the EFG. 

Can we understand the two slopes in Fig.^b)? Obvi- 
ously, n Pi should decrease with decreasing NN: Cd has no 
native 5p-electrons, hence in the free atom limit (NN=0) 
all n Pi should be zero. The reason why n Px (n Py ) grows 
faster with NN than n Pz can be understood from look- 
ing at the partial p x - and p 2 -D0S for two extreme cases 
(NN=3 and NN=9, see Fig-EJ- The p 2 -D0S has two dis- 
tinct contributions below Ep: a peak about 5 eV below 
Ep which nicely correlates with the s-DOS of Pd neigh- 
bors underneath the Cd atom, and hence indicates a ver- 
tical Pd-Cd s-p hybridization, and a second peak in the 
interval [—2, 0] eV that corresponds with the Pd d-DOS 
(Pd-Cd d-p hybridization). This second peak has almost 
disappeared for NN=3, which indicates that it originates 
mainly from Pd-neighbors in the horizontal plane. For 
the p x -DOS, there is only one peak, correlated mainly 
with the Pd d-DOS of neighbors in the horizontal plane. 
When going from NN=9 to NN=3, the number of Pd 
neighbors in the horizontal plane quickly drops (6 for 
NN=9, 4 for NN=8, 2 for NN=7 and for lower NN), 
while the number of Pd neighbors underneath Cd stays 
rather constant (3, 4 or 5, with the largest value for inter- 
mediate NN). Hence, n Px should quickly drop while n Pz 
tends to stay constant, which is exactly what we observed 
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FIG. 1: (a) Black symbols: calculated components of the EFG (see Tab. Hill black symbols refer to the right scale). White 
symbols: generalized asymmetry counts (see text), calculated from the n Vi in Tab. Hill (white symbols refer to the left scale). 
Gray symbols (for (110) surfaces only): average of An p and Arip, which would be the value of Arip—Arip in the hypothetical 
case of axial symmetry (77 = 0). Circles, triangles and squares are relative to the x,y,z components respectively of the plotted 
quantities. Dotted-dashed lines: An p obtained from the least squares fit in Fig.^b). The dotted lines are calculated using the 
extrapolation of Fig. 0^b). (b) The n Vi from Tab. Hill Gray symbols (for (110) surfaces only): average of n Px and n Py in the 
hypothetical case of axial symmetry (rj = 0). Dotted-dashed lines: least squares fit, with tentative extrapolation (dotted line) 
to NN=0 (free atom). Both in (a) and (b) the data for the xy-plane and for the z-direction are displayed separately for clarity. 
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FIG. 2: Partial p z - and p^-DOS for unrelaxed Cd on Pd 
(NN=3 and NN=9). 



in Fig. ^b). Reintroducing r\ ^ in our discussion only 
redistributes the charges within the horizontal plane, and 
therefore does not affect our conclusion. 

We emphasize that in contrast to the claimed sym- 
metry independence of the HFF (to be discussed in 
Sec. IIV If) , the coordination dependence of the p- 



distribution for sure depends on the exact spatial ar- 
rangements of the Pd neighbors. The behavior is so 
smooth only because we are dealing with fairly 'simi- 
lar' environments (Pd underneath, vacuum above). Fur- 
thermore, a tentative extrapolation of Fig. ^b) to very 
low coordination sites (NN=1, 2) suggests a negative V zz 
there. This is confirmed for the NN=2 environment for 
Cd on Ni to be discussed in Tab. IIIII where we calculate 
V zz = -2.4- 10 21 V/m 2 . 

In contrast to the discrete behavior suggested by the 
practical rule ('terrace site means large V zz , adatom site 
means low V zz ), we conclude that there is a continuous 
evolution as a function of the coordination number. The 
(100) and (111) terrace sites are both highly coordinated, 
and therefore have indeed the largest V zz , while the op- 
posite is true for the (100) and (111) adatom sites. The 
(110) terrace and adatom sites have intermediate coordi- 
nation - and therefore comparable V zz s - which is why 
the 'rule' does not seem to work for (110). However, there 
is nothing exceptional about the (110) surface. 



IV. MAGNETIC HYPERFINE FIELDS 

We now limit ourselves to Ni surfaces again and discuss 
the results obtained for the hypcrfine fields. Two main 
questions will be addressed. First, we want to verify 



and understand the experimentally proposed parabolic 
dependence of the Cd HFF on the number of Ni neigh- 
bors in the first coordination shell (NN), and we exam- 
ine whether and how the spatial arrangement of these 
neighbors influences the HFF. As a second problem, we 
broaden the scope to the entire 5sp series, placing 5sp 
impurity atoms only in NN=4 (adatom) and NN=8 (ter- 
race) coordinated sites at the Ni(100) surface. This will 
allow us to generalize the behavior observed for Cd to 
other probe atoms, and to propose qualitative explana- 
tions on the mechanisms ruling the observed HFF's. 



1. Cd HFF's on Ni surfaces 

In this section we focus on Cd probe atoms, placed 
in terrace and adatom positions at the three low-Millcr- 
index surface of fee Ni. This gives access to 6 differently 
coordinated sites, in addition to the fully coordinated 
substitutional bulk site and to a more artificial bridge 
site with NN=2. In order to test the sensibility of the 
HFF on the details of the cells, we performed also sev- 
eral calculations by varying the cell size and number of 
Ni layers in the slab. All the calculated Cd HFF for the 
differently coordinated systems are summarized in Ta- 
ble ]W\ Details on the cell are given in the format: (2D 
cell size , number of Ni layers). The Cd-Cd distance 
through the vacuum spacer is given in parentheses be- 
sides the HFF values. The long hyphen ' — ' is used to 
label cases that were not calculated. The NN=2 coordi- 
nated site has been achieved by placing the Cd atom in 
a non-crystallographic adatom position on Ni(lll). An 
'ideal unrelaxed position' has therefore no meaning for 
Cd on this site, which is indicated by the label 'n.p.' ('not 
possible'). For the fully coordinated bulk site, (NN=12) 
both a bulk and a slab calculation are reported, the lat- 
ter with the Cd placed in the middle layer of the Ni slab. 
The experimental results of Potzger— are reported in the 
last column. We first discuss the unrelaxed calculations, 
where all the atoms (including Cd) sit at their ideal fee 
position. The calculated values predict for the Cd HFF 
a change in sign for mid coordination and large positive 
values for low coordination, in agreement with the ex- 
perimental assignments. Changing the size of the slab, 
by adding Ni layers or by increasing the extension of the 
cells in the surface plane, results into some scattered val- 
ues which lay, except for NN=5 and to a lesser extent also 
for NN=4, within the aforementioned expected precision 
(± 5T). For the NN=5 case, i.e. adatom position at the 
Ni(llO) surface, the HFF's are found to attain large neg- 
ative values for some of the considered cells. We will 
come back to these puzzling results later on. 

As already discussed before, lattice relaxations are ex- 
pected in such open systems and might induce important 
changes to the HFF's. In fact it has been shown recently 
that their inclusion improves the agreement with the ex- 
perimental data, in the case of 5sp and 6sp impurities 
in bec bulk Fe>i2iii Due to the larger atomic volume of 
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FIG. 3: Calculated coordination number (NN) dependence of 
the Cd hyperfine field. Details are explained in the text. 



Cd with respect to Ni, an outward relaxation away from 
the surface is expected for both terrace and adatom po- 
sitions, for all low-index surfaces. Our calculations show 
Cd displacements from ideal fee positions towards the 
vacuum, that range from 0.60 to 0.90 A for the terrace 
atom and from 0.09 to 0.26 A for the adatom, depending 
on the surface orientation (Tab. 0. Displacements for 
the terrace site are larger than for the adatom site: the 
Cd atom strives for a Cd-Ni bond length of about 2.65 
A that is somewhat larger than the ideal fee Ni-Ni bond 
length of 2.48 A. Starting from an adatom position, this 
can be realized with less displacement. Minor relaxations 
appear also in the Ni atoms around the impurity. As ev- 
ident from Tab. IIVI the correction due to the relaxation 
is mostly within 7 T (except for NN=4 with corrections 
up to 35 T, see later). The correction is negative for 
the highest coordination numbers (NN=9 and 12) and 
positive for all others. 

For a better comparison, the theoretical results for 
the unrelaxed and relaxed systems are plotted in Fig. |3| 
together with the experimental results. When more than 
one value exists in Tab. IIVI the ones relative to the cell 
with the largest volume are selected, except for the more 
sensitive NN=4 and NN=5 cases where the values closer 
to the experiments are chosen (this will be justified be- 
low). The chosen HFF's are given in bold in Tab. IIVI 
As evident from Fig. |3 upon relaxation of the atoms in 
the cell the barycenter of the HFF curve moves towards 
more positive values for NN > 9, inducing a sign change 
already for NN = 7 (Cd on a Ni(110) terrace site). The 
largest changes are seen for NN=4 and NN=5, indicat- 
ing once more a high sensitivity of these environments. 
No clear overall improvement is found when relaxations 
are included, and the experimental data lays somehow 
between the relaxed and unrelaxed theoretical curves. 

In order to investigate the reason for the large varia- 
tions seen for the NN=4 and NN=5 coordinated sites, it 
is fruitful to look at the partial Density Of States (DOS) 



8 



TABLE IV: Hyperfine fields (T) of Cd in different Ni-environments. Experimental values are taken from Ref. ^] 











Hyperfine fields (in T) 






(J NN 


Type of cell 


non- relaxed 


relaxed 


Exp 


Bulk 


12 


(2x2x2) 


-9.8 


-12.0 


-6.9 


(100) Bulk 




(2x2, 7L) 


-10.0 


-11.7 




(Iff) Terrace 


9 


(75 x y/2 , 7L) 


-7.4 (12.18 A) 


-13.7 (10.49 A) 


-6.6 


(100) Terrace 


8 


(2x2, 5L) 


-9.8 (7.02 A) 


-1.9 (5.80 A) 


-3.5 






(2x2, 7L) 


-8.3 (7.02 A) 


-6.7 (5.82 A) 








(v^xV2, 7L) 


-3.6 (13.19 A) 


-5.9 (11.88 A) 




bulk-like ('random') 




(2x2x2) 


-9.9 






bulk-like (terrace) 




(2x2x2) 


-8.2 






bcc-bulk 




(2x2x2) 


-8.9 






(110) Terrace 


7 


(2 x 72 , 7L) 


-2.9 (7.45 A) 


1.1 (5.66 A) 








(2 x 272 , 9L) 


-6.8 (7.45 A) 







(110) Adatom 


5 


(2 x 72 , 5L) 


-4.8 (7.45 A) 


11.9 (6.91 A) 


4.3 






(2 x 72 , 7L) 


-14.2 (7.45 A) 


-6.5 (6.93 A) 








(2 x 272 , 7L) 


-21.7 (7.45 A) 







(100) Adatom 


4 


(2x2, 5L) 


-6.5 (7.02 A) 


-3.6 (6.51 A) 


7.3 






(2x2, 7L) 


3.4 (7.02 A) 


38.7 (6.63 A) 








(72 x 72 , 5L) 





15.9 (9.35 A) 








(72 x 72 , 7L) 


2.6 (9.67 A) 


22.9 (9.28 A) 








(72 x 72 , 9L) 


— 


31.9 (9.35 A) 








(72 x 72 , 11L) 




25.1 (9.35 A) 




bulk-like ('random') 




(2x2x2) 


-5.6 






bulk-like (adatom) 




(2x2x2) 


8.0 






bulk-like (free layer) 




(2x2x2) 


-4.2 






bee bulk like (adatom) 




(2x2x2) 


-9.6 






(111) Adatom 


3 


(272 x 272 , 5L) 


24.9 (8.11 A) 


31.9 (7.92 A) 


16.0 


(111) Adatom 


2 


(72 x 72 , 7L) 


n.p. 


33.4 (8.26 A) 





TABLE V: Cd perpendicular relaxations (in A) from ideal 
fee crystallographic site for the three low-index Ni surfaces. 
Everywhere Cd moves towards the vacuum. 



Site 


(100) 


(110) 


(111) 


Adatom 


0.20 


0.26 


0.09 


Terrace 


0.60 


0.90 


0.84 



of Cd with s-symmetry (only s-electrons contribute to the 
HFF). In Fig. 0] the majority and minority s-DOS of a 
Cd adatom on Ni(100) is shown in a small energy window 
in the vicinity of the Fermi energy. Since it is not clear 
a priori whether the problem arises due to the adopted 
approximation for the exchange-correlation functional, 
LDA-calculations are presented as well. The s-DOS 
shows several structures, among them a pronounced peak 
which for the majority channel lays right at Ep while it 
remains above Ep for the minority channel. Small vari- 
ations in the details of the cell (number of layers, 2D 



cell size, relaxation or not, . . .) or in the computational 
method (LDA/GGA, APW+lo or KKR, . . .) will push 
this peak in the majority channel below or above Ep. 
Since only the majority spin is involved, one has a net 
change of the s spin magnetic moment, and - because 
this is roughly proportional to the HFF - also a change 
of the HFF itself. For instance, in the two examples in 
Fig-El the s spin moment and HFF are 2.02- 10 -2 \xb and 
38.7 T for GGA, and 1.05 ■ 10~ 2 ^ B and 20.4 T for LDA. 
We pinpoint here that the LDA does not generally solve 
this sensitivity problem: in the NN=5 adatom unrelaxed 
case, for instance, similar negative s moments and large 
negative HFF are obtained both by LDA and GGA. As a 
general rule, whenever such a peak is observed so close to 
Ef, an enhanced sensitivity of the HFF is foreseen, which 
adds to the inner precision of the calculations. For NN=4 
and NN=5, we do observe such a peak in the s-DOS, for 
the other environments this is not the case. This explains 
the instability and wide scattering of the HFF's for those 
two coordination numbers in Tab. IIVI (note that this in- 
stability does not affect the calculations for V zz , which 
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FIG. 4: Partial s-DOS for relaxed Cd on Ni(100) calculated 
in the LDA and GGA approximations. 



arc ruled by p- rather than s-electrons). For impurity 
elements other than Cd, it might be that such sensitivity 
does not show up at all, or - if it does - it might do so for 
other values of the coordination number, depending on 
the details of the s-DOS close to Ep. The experimental 
values in Tab. IIVI are up to 20 T below the calculated 
values for NN>4. This indicates that in nature the ma- 
jority s-peak is above Ep, while in our calculations it is 
below. 

We now turn to the experimental parabolic NN- 
counting rule from Fig. [3J Both our relaxed and unre- 
laxed calculations show roughly the same trend (Fig. 01, 
and agree with the experimental trend. But how abso- 
lutely does this 'parabolic rule' holds? Is it really true - 
as concluded by Potzger et aliiS - that knowledge of the 
NN coordination is enough to predict the HFF on Cd in 
Ni-environments? In order to answer this question, we 
exploit the advantage of ab initio calculations that one 
is not restricted to environments that necessarily have 
to exist in nature. We can easily create artificial Cd- 
in-Ni-environments with an arbitrary number of nearest 
neighbors. That allows us to test the NN-counting rule in 
more situations than are experimentally accessible. The 
environments we created are 2x2x2 supercells for Cd in 
Ni, with a given amount of nearest neighbor Ni atoms 
removed (such that vacancies remain) . This removal was 
done pair by pair, with the requirement that the remain- 
ing cell still has inversion symmetry (out of the several 
possibilities for every NN, we calculated only one). This 
requirement makes such environments essentially differ- 
ent from the corresponding surface environments with 
the same NN, because the latter inevitably do not have 
inversion symmetry. In a first series, we start from the 
relaxed Cd-in-Ni supercell (labelled as 'semi-relaxed', be- 
cause after removal of the Ni atoms no further relaxation 
is done), in a second series all atoms are at the ideal 
Ni bulk positions. The results of both series are almost 
identical, and are reported in Fig. [3J Clearly, the same 



trend as in experiment and as in the surface calculations 
is present in the artificial bulk calculations. This proves 
that there is a basic truth in the NN-counting rule. On 
the other hand, the large difference between the bulk and 
surface calculations - especially for the low-coordination 
environments - is a clear sign that contributions from 
higher coordination shells are not negligible. 

If the requirement of inversion symmetry is removed, 
we can bridge the gap between these bulk-like cases and 
the surface slabs by calculating bulk-like cells where the 
first coordination shell is exactly the same as on a specific 
surface site. We did this for NN=8 and NN=4 fTab.Hv! 
the environment labelled as 'random' is the one with 
inversion symmetry). For NN=8, we could simulate in 
this way an environment that has exactly the same NN- 
coordination as the (100) terrace site. If the NN-counting 
rule were absolutely valid, we would find exactly the same 
HFF in the random and terrace-like bulk case. The re- 
sults are indeed quite close: -9.9 T and -8.2 T, which 
are values that are also not far from the slab calculation 
(-8.3 T). An even more daring test of the counting rule, 
is to put Cd at a substitutional site of an unrelaxed hypo- 
thetical bee Ni 2x2x2 supercell (16 atoms), with a lattice 
constant chosen such that the Ni-Ni distance (and hence 
also the unrelaxed Cd-Ni distance) is the same as in the 
fee case. Even in this very different kind of environment, 
the calculated HFF of -8.9 T follows the simple count- 
ing rule. All this is different for NN=4. Apart from the 
random environment (-5.6 T), we tested a configuration 
that is identical to the fee (100) adatom case (8.0 T), 
a configuration with all 4 Ni plus Cd in the same plane 
(a free layer, -4.2 T), and in a bee cell a configuration 
that is identical to a bee (100) adatom (-9.6 T). These 
4 numbers prove that - even in /cc-based environments 
only - the exact spatial configuration of the Ni neighbors 
in the first coordination shell can be important, leading 
to differences of more than 13 T. It is not surprising to 
find this effect for NN=4 rather than NN=8, the for- 
mer being identified before as a sensitive case. From this 
analysis we conclude that although the NN-counting rule 
certainly indicates a trend, there can be substantial devi- 
ations from it for specific environments. In such sensitive 
environments, the spatial arrangement of the neighbors is 
important as well. There is probably some luck involved 
that for Cd in Ni the behavior in nature is so smooth 
as experimentally observed, and there is no fundamental 
reason why the data could not have been considerably 
more scattered around the parabolic trend. 



2. HFF's of the 5sp series 

As a last part of this study we now present a survey for 
the 5sp impurities from Cd to Ba in Ni-environments, to 
see if and how the NN-counting rule can be extended 
to other impurities. We use the strategy applied by 
Mavropoulos et for 4sp impurities on Fc and Ni 

surfaces, and very recently and independently from this 
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FIG. 5: Hyperfine fields of the whole 5sp series (Cd -^Ba 
in bulk Ni (circle, NN=12), at a terrace position (squares, 
NN=8) and at an adatom position (diamonds, NN=4) on the 
Ni(100) surface. 




FIG. 6: Cartoons for the majority and minority partial s-DOS 
of (a) a 5sp impurity in bulk Ni and (b) a 5sp impurity at a 
Ni surface. The vertical line indicates the Fermi energy, the 
name of the elements indicates for which element a particular 
picture is representative. This picture is inspired by Ref. 0. 



work also for 5sp impurities on Ni surfaces wL We take the 
bulk environment (NN=12) and the terrace (NN=8) and 
adatom (NN=4) environments for the (100) surface, and 
calculate the HFF for the 9 elements from Cd to Ba in 
those environments. The cell sizes chosen for those cal- 
culations are the ones labelled as in Tab. II VI as (2x2x2), 
(2x2, 5L) and (2x2, 7L), respectively. Because we 
are looking again for gross trends now and not for fine 
details, relaxations were not included. The results for 
the 5sp HFF's are displayed in Fig. |3J Exactly the same 
behavior as Mavropoulos et al. observed for 4sp and 5sp 
impurities is seen here, which mutually supports the va- 
lidity of the very different computational methods that 
were used. In the bulk environment, the HFF starts at 
about -10 T for Cd, strongly increases with increasing 
atomic number Z, and reaches a maximum near the mid- 
dle of the series: ~46 T for I. Then it decreases again, 
and at the end of the series turns back to values close 
to -10 T. When the coordination number is reduced, the 
main peak of the HFF curve moves to heavier elements 
and an additional structure - that for NN=8 is more a 
broad shoulder than a peak - appears at the beginning 
of the series. For NN=4 two clear structures are evident, 
and the HFF increases and decreases (with less intense 
variations than for the bulk) twice in the course of the 
5sp series. 

The microscopic origin of Fig.JSJcan be understood by 
a slight extension of arguments given by Mavropoulos et 
a liiii for 4sp impurities. Later on, we will then derive 
(and test) generalized NN-counting rules for all 5sp im- 
purities in Ni environments from them. For a systematic 
explanation, let us go back to the origin of hyperfine fields 
in ferromagnets (see Ref. for a detailed and instructive 
review). As the hyperfine field in our cases is dominated 
by the Fermi contact contribution, we have to care about 
the details of the bond between the 5s states and its en- 
vironment (here Ni-3d). It has been known for a long 



tim o 54 i 55 that in the case of such an s-d bond the local 
s-DOS of the impurity shows a characteristic depression 
a few eV below the fermi energy: the 'antiresonance dip' 
(AR). The position of the AR is mainly determined by 
the host material (Ni), and not by the impurity. The 
states below the AR arc bonding states, the states above 
are antibonding. The up and down states are exchange 
split, such that at first sight one expects an excess of 
s-up over s-down, resulting in a positive hyperfine field. 
Due to a different s-d hybridization for up and down elec- 
trons, however, the number of s-up below AR will be di- 
minished, while the number of s-down will be enhanced. 
Above AR, the situation is opposite^ The final result is 
that the impurity s-moment (and hyperfine field) will be 
negative in the beginning of the sp-series, where the ef- 
fect of the bonding states is dominant (Fig.^jJ-a-l). In the 
second half of the series, also the antibonding states will 
get filled, and because they have to be squeezed between 
AR and the fermi energy, they have to develop a sharp 
peak in the DOS. The exchange splitting of this peak is 
responsible for the large positive HFF at the end of the 
sp-series (Fig.|SJ-a-2/3), which quickly drops to small and 
negative values again if also the down antibonding states 
are below the Fermi energy (Fig.|J3-a-4). Mavropoulos et 
al. have shown by group theoretical arguments that in 
the case of reduced point group symmetry for the impu- 
rity (as on surfaces), the antibonding part of the impu- 
rity s-DOS is split in two parts. This splitting is more 
pronounced if the impurity is in a more non-bulk-likc en- 
vironments, i.e. it is more pronounced for NN=4 then for 
NN=8. Let us take the NN=4 case with a clear splitting. 
When going from Cd to Ba, we evolve through the dif- 
ferent stages of Fig. EJ-b, which explains the double-peak 
structure of the HFF for NN=4 in Fig. 

An aspect of Fig. [5] that has not been discussed by 
Mavropoulos et al., is the physical origin of the coordi- 
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FIG. 7: (a) The first of the two antibonding s-peaks for ma- 
jority and minority spin, (b) The s-moment derived from (a) 
by subtracting - at a particular energy - the integral of the 
minority s-DOS up to that energy from the integral of the ma- 
jority s-DOS up to that energy. In both (a) and (b) 3 typical 
cases are drawn: high coordination = broad band width (thin 
full line), medium coordination = medium band width (thick 
full line) and low coordination = small band width (dashed 
line). The horizontal arrows indicate the region where the 
Fermi energy falls for the indicated elements (first half of the 
5sp-series). This story is can be repeated with the second 
of the two antibonding peaks, starting from I for which the 
Fermi energy falls at the place indicated by the vertical arrow. 



nation number dependence of the HFF for a particular 
element: why is e.g. for In the NN=4 HFF the larger one 
and the NN=12 the smaller one, while this is reversed for 
e.g. Te? This we will explain by the cartoon in Fig. 
The upper part of Fig. schematically shows the first 
of the two antibonding s-peaks of Fig. [5Jb, for any par- 
ticular impurity. If one lowers the coordination number 
of the impurity, the band-width of these peaks will de- 
crease - an obvious fact, which we clearly observe in our 
calculations. Fig.0-a shows the same situation for 3 typ- 
ical band-widths: large (NN=12), medium (NN=8) and 
small (NN=4). The bottom part of Fig. shows the s 
spin moment derived from Fig. 0a as a function of energy 
(found by subtracting the integral of the down-peak from 
the integral of the up-peak, where the integrals are made 
up to the energy under consideration). Everything now 
depends on where the Fermi energy lies in Fig.0 If it falls 
in the region indicated by Cd-In-Sn, Fig. 0a corresponds 
to Fig. [SJ-b-l. At the corresponding energy in Fig. 0b, 



the s-moment (HFF) for the small band-width (NN=4) is 
larger then for the medium band- width (NN=8), which 
in turn is larger than for the large band- width (NN=12). 
If the Fermi energy falls in the region indicated by Sb- 
Te, Fig. 0a corresponds to Fig. and the sequence 

of s-moments is reversed. After Te, this story repeats 
for the second antibonding peak (for I, the Fermi energy 
will be at the position marked by the vertical arrow, but 
of course in the second series of peaks), but will be in- 
creasingly less clear due to the presence of the 6s states 
that start to manifest themselves around the Fermi en- 
ergy. That is the reason for the more chaotic evolution 
for Cs and Ba. Of course, Fig. is a cartoon only, and 
its conclusions should not be taken too literally: the real 
DOS are not Gaussians as used in the cartoon, and there- 
fore the details of the hyperfine field evolution might be 
different. Nevertheless, it captures the basic mechanism. 
Summarizing, we conclude that the physical mechanism 
behind Fig. can be understood from a combination of 
three basic features: i) the double peak structure of the 
antibonding peaks, ii) the decrease in the band width 
and hence the increase of peak height - upon reduction 
of the coordination number, and iii) the position of the 
Fermi energy with respect to the peaks. 



We now take Fig.0as a source of inspiration to extend 
the parabolic NN-counting rule proposed by Potzger et 
al. (Ref. Il9l) for 5sp impurities other than Cd. It can be 
seen from Fig. that the Cd-HFF for bulk and NN=8 is 
almost the same and negative, while the value for NN=4 
is small and positive: this is the parabolic behavior seen 
in experiment. In the same way, we can then deduce that 
for In and Sn as impurities, the HFF should monotoni- 
cally rise from NN=12 to NN=4 (it is interesting to note 
that for the experimentally 'easily' accessible Mossbauer 
probe 119 Sn, Fig. suggests a linear behavior). Between 
Sn and Sb all lines cross, such that for Sb to I the HFF 
monotonically decreases from NN=12 to NN=4. For Xe 
to Ba, there is non-monotonic behavior instead. We have 
checked this deduction by calculating the artificial bulk- 
like (non-relaxed) environments as discussed before, but 
now for Te and Ba as impurities (both are taken as a rep- 
resentative for the region of monotonic decrease and the 
non-monotonic region, just as Cd is a representative of 
the region of monotonic increase). In these bulk- like cells, 
we can more easily create environments with NN differ- 
ent from 4, 8 and 12. The results are given in Fig. 
The general trends shown in Fig. are the same as the 
ones which could be inferred from the surface calcula- 
tions of Fig. [3J monotonic (parabolic) increase for Cd, 
monotonic decrease for Te and non-monotonic behavior 
for Ba. In the light of these results we therefore conclude 
that each of the 5sp impurities in Ni has its own typical 
coordination number counting rule. 
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FIG. 8: Coordination dependence of the HFF for selected, i.e. 
Cd,Te and Ba, 5sp elements as obtained by bulk calculations 
(see text for details). 



V. CONCLUSIONS 

We have undertaken an extensive confrontation be- 
tween state-of-the-art ab initio calculations and a high- 
quality data set - experimentally collected during the 
past 20 years - of electric-field gradients and magnetic hy- 
perfine fields of Cd at magnetic and non-magnetic metal- 
lic fee surfaces. We conclude that the experimental prac- 
tice for (100) and (111) surfaces of assigning a large V zz 
to terrace sites and a low V zz to adatom sites is cor- 
rect. However, this 'discrete' rule is nothing more than a 
manifestation of a continuous evolution of the EFG com- 



ponents as a function of coordination number, an evolu- 
tion for which we have pointed out the physical mecha- 
nism. With this insight, the behavior of the (110) surface 
is not exceptional at all. The experimentally suggested 
parabolic-like coordination number dependence for the 
HFF of Cd at Ni surfaces is confirmed as being a reli- 
able trend, but we warn that it is just a trend and not 
a rigorous rule: sensitive environments exist, for which 
the spatial arrangement of the Ni neighbors considerably 
influences the HFF. We have explained in detail the phys- 
ical mechanism behind the HFF for all 5sp impurities at 
Ni surfaces, by combining knowledge from the literature 
and new insight. In particular we have generalized the 
parabolic NN-counting rule for Cd to other 5sp impuri- 
ties, showing that each impurity has its own typical rule, 
and explaining why this is so. We hope to have demon- 
strated that ab initio calculations can greatly enhance 
the physical insight in an experimentally complex prob- 
lem. 

Acknowledgments 

One of the authors (V.B.) acknowledges fruitful discus- 
sion with Dr. K. Potzger. Part of the calculations were 
performed on computer facilities granted by an INFM 
project Iniziativa Trasversale Calcolo Parallelo at the 
CINECA supcrcomputing center. Another part of com- 
putational work was performed on computers in Lcu- 
ven, in the frame of project G. 0239. 03 of the Fonds voor 
Wetenschappelijk Onderzoek - Vlaanderen (FWO). The 
authors are indebted to L. Verwilst and J. Knudts for 
building and maintaining the pc-cluster in Leuven. M.C. 
thanks the Onderzoeksfonds of the K.U. Leuven for grant- 
ing a post-doctoral fellowship (F/02/010). Knowledge, 
help and advice from the active WIEN2k user commu- 
nity are warmly acknowledged. 



* Electronic address: Stefaan.Cottenier@fys.kuleuven.ac.be 
t Electronic address: Bellini.Valerio@unimore.it 

1 R. Vianden, in Nuclear Physics Application on Materials 
Science, edited by E. Recknagel and J. C. Soares (Kluwer 
Academic Publishers, 1988), p. 239. 

2 G. Schatz, in Nuclear Physics Application on Materials 
Science, edited by E. Recknagel and J. C. Soares (Kluwer 
Academic Publishers, 1988), p. 297. 

3 G. Krausch, R. Fink, K. Jacobs, U. Kohl, J. Lohnmuller, 
B. Luckscheiter, R. Platzer, B.-U. Runge, U. Wohrmann, 
and G. Schatz, Hyperfine Interactions 78, 261 (1993). 

4 H. Haas, Z. Naturforsch. 50a, 407 (1994). 

5 J. Voigt, R. Fink, G. Krausch, B. Luckscheiter, R. Platzer, 
U. Wohrmann, X. L. Ding, and G. Schatz, Phys. Rev. 
Letters 64, 2202 (1990). 

6 H. H. Bertschat, H.-H. Blaschek, H. Granzer, K. Potzger, 
S. Seeger, and W.-D. Zeitz, Phys. Rev. Letters 80, 2721 
(1998). 

7 H. Akai, M. Akai, S. Bliigel, R. Zeller, and P. H. Dederichs, 
J. Magn. Magn. Materials 45, 291 (1984). 



8 S. Bliigel, H. Akai, R. Zeller, and P. H. Dederichs, Phys. 
Rev. B 35, 3271 (1987). 

9 H. Akai, M. Akai, S. Bliigel, B. Drittler, H. Ebert, K. Ter- 
akura, R. Zeller, and P. H. Dederichs, Progr. Theor. Phys. 
Suppl. 101, 11 (1990). 

10 T. Korhonen, A. Settels, N. Papanikolaou, R. Zeller, and 
P. H. Dederichs, Phys. Rev. B 62, 452 (2000). 

11 S. Cottenier and H. Haas, Phys. Rev. B 62, 461 (2000). 

12 H. Haas, Hyperfine Interactions 129, 493 (2000). 

13 V. Bellini, R. Zeller, and P. Dederichs, Phys. Rev. B 64, 
144427 (2001). 

14 C. O. Rodriguez, M. V. Ganduglia-Pirovano, E. L. Peltzer 
y Blanca, M. Petersen, and P. Novak, Phys. Rev. B 63, 
184413 (2001). 

15 C. O. Rodriguez, M. V. Ganduglia-Pirovano, E. L. Peltzer 
y Blanca, and M. Petersen, Phys. Rev. B 64, 144419 
(2001). 

16 R. N. Nogueira and H. M. Petrilli, Phys. Rev. B 60, 4120 
(1999). 

17 S. Frota-Pessoa and S. Legoas, Hyperfine Interactions 133, 



13 



207 (2001). 

18 J. A. Gomez and D. Guenzenburger, J. Phys.: Condensed 
Matter 14, 12311 (2002). 

19 K. Potzger, A. Weber, H. Bertschat, W.-D.Zeitz, and 
M. Dietrich, Phys. Rev. Letters 88, 247201 (2002). 

20 M. J. Prandolini, Y. Manzhur, A. Weber, K. Potzger, 
H. H. Bertschat, and M. Dietrich, abstract presented at 
ICM2003. 

21 H. Haas, Hyperfine Interactions 84, 69 (1994). 

22 B. Lindgren, Hyperfine Interactions 34, 217 (1987). 

23 B. Lindgren, Europhys. Lett. 11, 555 (1990). 

24 B. Lindgren and A. Ghandour, Hyperfine Interactions 78, 
291 (1993). 

25 B. Lindgren, Z. Naturforsch. 57a, 544 (2002). 

26 Ph. Mavropoulos, N. Stefanou, B. Nonas, R. Zeller, and 
P. H. Dederichs, Phys. Rev. Letters 81, 1505 (1998). 

27 Ph. Mavropoulos, J. Phys.: Condensed Matter 15, 8115 
(2003). 

28 P. Hohenberg and W. Kohn, Phys. Rev. 136, 864 (1964). 

29 W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

30 P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and 
J. Luitz, WIEN2k, An Augmented Plane Wave + Local 
Orbitals Program for Calculating Crystal Properties, (Karl- 
heinz Schwarz, Techn. Universitat Wien, Austria) (1999), 
ISBN 3-9501031-1-2. 

31 J. P. Perdew, S. Burke, and M. Ernzerhof, Phys. Rev. Let- 
ters 77, 3865 (1996). 

32 A. Settels, T. Korhonen, N. Papanikolaou, R. Zeller, and 
P. H. Dederichs, Phys. Rev. Letters 83, 4369 (1999). 

33 H. M. Petrilli, P. E. Blochl, P. Blaha, and K. Schwarz, 
Phys. Rev. B 57, 14690 (1998). 

34 G. Kresse and J. Furthmiiller, Comp. Mater. Science 6, 15 
(1996). 

35 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

36 D. A. Shirley, S. S. Rosenblum, and E. Matthias, Phys. 
Rev. 170, 363 (1968). 

37 E. Fermi, Z. Physik 60, 320 (1930). 



38 P. Dufek, P. Blaha, and K. Schwarz, Phys. Rev. Letters 
75, 3545 (1995). 

39 T. Klas, J. Voigt, W. Keppner, R. Wesche, and G. Schatz, 
Phys. Rev. Letters 57, 1068 (1986). 

40 E. Hunger and H. Haas, Surf. Sci. 234, 273 (1990). 

41 J. Voigt, Ph.D. thesis, Universitat Konstanz (1990). 

42 R. Fink, B.-U. Runge, K. Jacobs, G. Krausch, 
J. Lohmiiller, B. Luckscheiter, U. Wohrmann, and 
G. Schatz, J. Phys.: Condensed Matter 5, 3837 (1993). 

43 R. Wesche, R. Fink, T. Klas, G. Krausch, R. Platzer, 
J. Voigt, and G. Schatz, J. Phys.: Condensed Matter 1, 
7407 (1989). 

44 R. Fink, R. Wesche, T. Klas, G. Krausch, R. Platzer, 
J. Voigt, U. Whrmann, and G. Schatz, Surf. Sci. 225, 331 
(1990). 

45 T. Klas, R. Fink, G. Krausch, R. Platzer, J. Voigt, 
R. Wesche, and G. Schatz, Surf. Sci. 216, 270 (1989). 

46 G. Fillebock, private communication to B. Lindgren. 

47 T. Klas, R. Fink, G. Krausch, R. Platzer, J. Voigt, 
R. Wesche, and G. Schatz, Europhys. Lett. 7, 151 (1988). 

48 E. Hunger, H. Haas, and H. Grawe, Hyperfine Interactions 
60, 999 (1990). 

49 P. Blaha, K. Schwarz, and P. H. Dederichs, Phys. Rev. B 
37, 2792 (1988). 

50 C. Ambrosch-Draxl, P. Blaha, and K. Schwarz, J. Phys.: 
Condensed Matter 1, 4491 (1989). 

jl K. Schwarz, C. Ambrosch-Draxl, and P. Blaha, Phys. Rev. 
B 42, 2051 (1990). 

52 S. Lany, P. Blaha, J. Hamann, V. Ostheimer, H. Wolf, and 
T. Wichert, Phys. Rev. B 62, R2259 (2000). 

53 S. J. Asadabadi, S. Cottenier, H. Akbarzadeh, R. Saki, and 
M. Rots, Phys. Rev. B 66, 195103 (2002). 

54 K. Terakura, Progr. Theor. Phys. 46, 1007 (1971). 

55 K. Terakura, J. Phys. F: Met. Phys. 7, 1773 (1977). 

56 E. Daniel and J. Friedel, J. Phys. Chem. Solids 24, 1601 
(1963). 



